Exact Chiral Spin Liquid with Stable Spin Fermi Surface on the Kagome Lattice 
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We study an exactly solvable quantum spin model of Kitaev type on the kagome lattice. We find 
a rich phase diagram which includes a topological (gapped) chiral spin liquid with gapless chiral 
edge states, and a gapless chiral spin liquid phase with a spin Fermi surface. The ground state of the 
current model contains an odd number of electrons per unit cell which qualitatively distinguishes it 
from previously studied exactly solvable models with a spin Fermi surface. Moreover, we show that 
the spin Fermi surface is stable against weak perturbations. 
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Introduction-The search for exotic ground states of 
quantum many-body systems has been at the fore in 
experimental and theoretical research for a number of 
decades, yet nature continues to tantalize us with ever 
more interesting behaviors. In recent years, systems with 
some type of topological order beyond Landau's broken 
symmetry principle [TJ [2] , such as the fractional quantum 
Hall effect [3] and quantum spin liquids |3] , have garnered 
a great deal of attention. A quantum spin liquid (QSL) 
is an insulating state that does not exhibit any conven- 
tional magnetic order at zero temperature, typically due 
to a delicate balance of competing and/or "frustrating" 
effects (for a review, see Ref. [HI E] ) - 

Various flavors of QSLs with different types of quasi- 
particle excitations and braiding statistics have been pro- 
posed [JHS]- However, the reliability of the approxima- 
tions typically invoked in their solution are often im- 
perfectly understood. Thus, it is useful to find exactly 
solvable (even though often contrived) models with spin 
liquid ground states so as to establish their stability as 
phases of matter and to understand their essential phys- 
ical properties. The existence of such states is com- 
pellingly established [10] for quantum dimer models [IT] 
on the triangular lattice. 

More recently, the discovery by Kitaev [T2] of an 
exactly solvable, interacting two-dimensional spin-1/2 
model on the honeycomb lattice with a spin liquid phase 
came as a great boon to the study of exotic ground states 
in magnetic systems. Since then, variants [13 of Kitaev's 
model on trivalent lattices have appeared which realize 
chiral spin liquids with non-Abelian anyons [T3] and spin 
liquids with a spin Fermi surface [T5HT7] . They have also 
been used to study quantum critical points [HI [19] , en- 
tanglement entropy and entanglement spectrum [20] , and 
edge solitons [2JJ among others [22J 123] • If the coordina- 
tion number of the lattice is larger than three, a general- 
ization of the Kitaev model, the so-called T-matrix model 
[T51 I24T[2"5] , with extra degrees of freedom on each site is 
needed to ensure exact solvability. These additional de- 
grees of freedom can be either interpreted as an "orbital" 
degree of freedom [2SI 123 , or one may consider each site 



as having a larger spin, but no orbital degrees of freedom 
01 [30]. 

Partly motivated by the recent interest in finding a spin 
liquid with a spinon Fermi surface in Herbertsmithite-a 
Mott insulator on the kagome lattice [31] -we study an ex- 
actly solvable spin-3/2 model on the kagome lattice (see 
Fig. [I]). Our main finding is that our model possesses a 
gapless spin liquid phase with a finite spin Fermi surface. 
We show that the spin Fermi surface is stable against 
any weak perturbations. The stability follows from the 
lack of inversion symmetry, 7r rotation and time-reversal 
symmetry (that is spontaneously broken) in the ground 
state. Both the odd number of electrons per unit cell in 
the ground state and the stability of spin Fermi surface 
qualitatively distinguishes our model from other known 
exactly solvable models with a finite spin Fermi surface 
[T5HT7] . We also find a gapped Abelian chiral spin liquid 
phase [TJ], which possesses two gapless chiral Majorana 
edge states for a lattice geometry with boundary. 

Model- The T-matrix model describing 5 = 3/2 spins 
[T5] we study respects the translational and threefold ro- 
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FIG. 1. (a) is a schematic representation of kagome lat- 
tice which is a network of corner sharing triangles in two- 
dimensions with 3 sites per unit cell, (b-e) are 4 energetically 
distinct flux configurations preserving the translational sym- 
metry of the lattice. A doubled (magnetic) unit cell is required 
for (d) and (e) where the total flux through a unit cell is ±7r. 
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tational symmetry of the kagome lattice and is given by 

n = J ^ E r ! r ? + ^v £ rfr 4 + J 5 £r 4 5 
+j A £ rfrf + £ rfrf, (i) 



(a) 



where we have distinguished the nearest neighbor cou- 
plings Jij as Ja, Jv and as J A and Jy based on 
whether the link (ij) belongs to an up (A) triangle or 
down (V) triangle, and (ij) is taken in the counter- 
clockwise sense for each triangle. Locally, the five F- 
matricies satisfy a Clifford algebra, {Tf,T b } = 25 ab , 
where a,b = 1,...,5, and T ab = [T a , T b ]/(2i). In terms 
of the components of the spin S — 3/2 operators [T5I I30j . 



-^={S Z ,S X },T 3 
V3 



r 4 = ~[(s*) 2 -(^) 2 ], r 5 = (^) 2 -^ 



{S X ,S V }, 



(2) 



With the identification ([2]), it is clear the model ([!]) 
has a global Ising spin symmetry under 180° rotations 
about the z-axis, and possesses time-reversal symmetry 
(TRS), although TRS will be spontaneously broken in 
the ground state as we describe below, in addition to the 
translational and threefold rotational lattice symmetry 
mentioned above. 

Method of Solution-The key to the exact solvability 
02] of the model ([!]) is an infinite number of conserved 
operators, W p , that satisfy [Tl^,^] = and [Wp,Wp/' 
for all elementary plaquettes labeled by p and p' . As 
shown in Fig. [l] the kagome lattice contains elementary 
up and down triangular plaquettes, and hexagonal ones 
Specifically, 



rfrfrf 4 , and 



W = r 23 r} 4 r 23 r ; 14 r^ri 4 , where the sites are taken in 
a counter-clockwise fashion and for Wq the first link (ij) 
is assumed to lie on a A. 

In order to solve the Hamiltonian we introduce a 
Majorana representation of the T-matricies [13] , 



a5 



(3) 



with a = 1,2,3,4. There are thus 6 Majorana species 
on each site i: , £ 2 , £ 4 , Cj, c£j}. The Majorana rep- 
resentation enlarges the spin-3/2 Hilbert space, so that 
one must enforce the constraint Di — — r|r?r 3 r 4 r| = 

'i/s 2 i'; ! i,' = 1) namely for any physical state |4 , ) p hys 
phys — |*) P hys for any i. From the enlarged Hilbert 
space, physical states are obtained by applying the pro- 
jection operator: P = \[ i [ 1+ 2 Pi ] , where the product is 
over all sites in the lattice. 

Using the relations (J3|, the Hamiltonian becomes 



E 

(ij) 



[j l:j iUijCiCj + J'ijiuijdidj] + J 5 >J icidi, (4) 
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Fig 1(b) 
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Fig 1(e) 
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FIG. 2. (color online) (a) Ground-state energy per site of 
as a function of J5/J with Ja = Ja = Jv = Jv — J 
for the uniform flux configurations shown in Fig[TJb-e). (b) 
Ground-state flux configuration among Fig[TJb-e) as a func- 
tion of Ja , Jv , J5 for J A + Jv + Ji — 1 with Jv = Jv and 
Ja = Ja- The coloring scheme follows that of the legend in 
Figga). 



with Uij = — if ij G A and mj = 
ifoj € V. The original Hamiltonian is obtained 



where ity = 



through "H = PHP. Since [u„-,H] = [ 



} = 0, the 



original spin Hamiltonian has been reduced to a model 
of free Majorana fermions moving in a static background 
of Z2 gauge fields Q2]: in @ the can be replaced by 
their eigenvalues ±1. It remains to determine the pat- 
tern of the (up to gauge transformations) that yields 
the lowest energy state. We will limit our discussion to 
positive J couplings only. 

In terms of static Z 2 gauge fields, we define a flux <fi p 
] = via exp(i<j) p ) = Yiikep^jk) wnere jk is taken counter- 
clockwise, on each elementary plaquette p. It is clear 
that 4> p = ±7r/2 in the triangular plaquettes, <f> p — 0,7r 
in the hexagonal plaquettes. Since W p = T[(ij)e P u iji 
where ij is also taken counterclockwise, exp(i(f> p ) = — W p 
(— iWp) when p is hexagon (triangular) plaquettes. Un- 
der time reversal symmetry, W p — > ±W p where +(— ) is 



for hexagon (triangle) plaquettes; it follows that <p p — ¥ 
—tfip for triangle plaquettes while 4>p remains unchanged 
for hexagon plaquettes. Consequently, a ground state 
with a certain flux pattern {4> p } (regardless of any par- 
ticular choice of ) spontaneously breaks time reversal 
symmetry and its energy must be degenerate with the 
ground state with flux pattern obtained from {4> p } by 
changing (j) p — > — <p p on all triangular plaquettes |14j . 

In general, the nature of a ground state and its exci- 
tations crucially depends on the flux pattern {4> p } in the 
ground state. As we will see below, the richness of this 
model is due in part to the freedom to add additional 
terms to that preserve the exact solvability, but can 
be used to select different (gauge-inequivalent) configu- 
rations of the ground state. 

Chiral Spin Liquid-A chiral spin liquid (CSL) breaks 
time-reversal symmetry spontaneously and has no con- 
ventional magnetic order JS] Q31 [35] ■ As we have dis- 
cussed, the current model (JIJ indeed breaks TRS spon- 
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FIG. 3. (a) Band structure on a cylindrical geometry for Ja = 
J' A — 1.0, Jv = Jv = 0.8, J5 = 0. There are two gapless 
chiral Majorana edge states (dotted lines) which overlap on 
each other, (b) For J A = 1.0, J' A = 0.6, J v = 0.9, Jv = 0.5, 
and J5 = 0.1, the two gapless edges states separate. These 
ground states are thus CSLs with a spectrum Chern number 
(±2) and the vortices obey Abelian statistics. 



taneously. It remains to be shown that a ground state 
exists without breaking spin-rotational or translational 
symmetry. To determine the ground-state flux configu- 
ration, we compute the ground-state energy Eq of the 
Hamiltonian Q as a function of the Uy. 

For a bipartite lattice, the ground state flux pattern 
is determined by the Lieb theorem |33j and the flux is 
always uniform. However, the kagome lattice is not bi- 
partite, so we must determine it numerically and ana- 
lytically in some special limits [2]. By diagonalizing 
systems of up to about 10 4 sites and computing Eq, 
we indeed find that for various possible parameters in 
Hamiltonian (1), the ground state has a uniform flux pat- 
tern. There are four possible time-reversal-inequivalent 
uniform flux configuration labeled by {4>a, 4>v, 0o} = 
{f,f,7r},{-f,f,7r},{f,f,0},{-f,f,0} as shown in 
Fig. l(b-e). The ground state flux pattern depends 
on the relative size of all the coupling constants in 
Eq.([T]), as seen in Fig. [2j For instance, when J 5 <C 
{Ja, Jy , J' A , J^}, a uniform flux state of {f,f,0} or 
{— f j — f ; 0} is favored. 

For small J 5 , the ground state has the uniform flux con- 
figuration {-, |,0} shown in FigJTJd) and all its excita- 
tions are fully gapped for a system with periodic bound- 
ary conditions. The total flux per lattice unit cell is 7r; 
it follows that the magnetic unit cell must be doubled by 
a gauge choice. Note that even though a specific gauge 
choice Uij with this flux configuration would "seemingly" 
double the unit cell and break the threefold rotational 
symmetry of the lattice, the physical ground state ob- 
tained after the projection is translationally and three- 
fold rotationally invariant because the gauge choice obeys 
the corresponding projective symmetry group transfor- 
mations that combines a physical symmetry (here, trans- 
lation or threefold rotation) and an appropriate gauge 



FIG. 4. The Fermi surface (solid line) for the 
flux configuration {f,f,vr} and { Ja, Jv, J'a, Jv; J$} = 
{1.0,0.3,0.8,0.5,1.4}. The dashed hexagon is the Brillouin 
zone boundary. Note that there is only one Fermi pocket for 
this set of parameters and the three pockets shown are related 
by a reciprocal lattice vector and are thus equivalent. 

transformation 34 . Consequently, no symmetries other 
than the time reversal arc broken in the physical ground 
state. 

For the system with cylindrical symmetry, Fig. [3^a) 
and (b) illustrate the existence of chiral edge modes 
indicating a finite spectral Chern number of ±2. For 
Ja. = J' A , Jv = J\" anc ^ ^5 = 0) t ne c and d Majorana 
fermions decouple and have identical spectrum. It follows 
that the two gapless chiral Majorana edge states overlap 
as shown in Fig. [3] (a). Generically, the two chiral edge 
modes are separated, as shown in Fig. [3^b) . Moreover, 
the arguments in Refs. [T5J [23 [35] for spin-correlations 
apply here as well-the spin correlations vanish beyond 
nearest neighbors, establishing a CSL phase with a finite 
spectral Chern number and gapless chiral edge states on 
the boundary as a ground state of ([I]). Because of the 
even spectral Chern number (±2), we expect the vortex 
excitations are Abelian. Nonetheless, if the parameters 
in the Hamiltonian doubles the unit cell, a gapped non- 
Abelian chiral spin liquid can also be obtained [36 . 

Gapless spin liquid with stable spin Fermi surface- For 
arbitrary values of Ja , Jy , J'a-i J\j in > we found that 
the system is always gapped so long as we kept J5 = 0, 
where W p = — 1 for all p |36j . By increasing J5 , we found 
that the system closes its gap at some critical values of 
J5 and the ground state flux configuration changes when 
a critical point is crossed [see Fig. [2f[a)]. Note that the 
energy differences between different flux states are sur- 
prising small, as shown in Fig. [2ja). For a generic J5 not 
at those critical points, the system stays gapped, which 
seems to make the search for a spin Fermi surface difficult 
for the Hamiltonian (1). Fortunately, while changing J5 
a particular flux configuration can be always favored by 
adding a pure flux term that does not destroy the exact 
solvability of the model. Specifically, to stabilize a flux 
configuration we add the following term to Eq. 

n FS = -c^VFo-/3 J2 W A Wv (5) 

O (A,V) 

The required value of a and (3 to stabilize a desired flux 
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pattern depends on J 5 , but is typically rather small [THl 
[17] as the energy differences between different flux states 
are small-see Fig. [2ja). 

For appropriate positive a and /3, the ground state 
with the flux configuration {§, §,7r} (or {— |,— §,7r}) is 
favored. For this flux configuration, we obtain a finite 
Fermi surface as shown in Figg (a) for J A = 1.0, J v = 
0.3, J' A = 0.8, = 0.5, and J 5 = 1.4. The dashed line 
of the hexagon is the Brillouin zone boundary. Note that 
there is only one Fermi pocket around one inequivalent 
corner for this set of parameters. The Fermi surface has 
a C3 rotational symmetry around the zone corner that is 
expected from the C3 symmetry of the model. Varying 
the parameters in Eq. (1) by a small amount will only 
change the shape and size of the Fermi surface but the 
Fermi surface itself is robust. 

An important question to ask about a Fermi surface is 
whether it is stable against any weak perturbations. In 
2D, a generic Fermi surface can not be fully gapped by 
a weak "density wave" due to the lack of perfect nest- 
ing. However, a generic Fermi surface can still be fully 
gapped by a "pairing" term with momentum Q if the 
spectrum at k and — k + Q are degenerate for any k and 
some fixed Q. In general, time reversal, tt rotation, or 
lattice inversion symmetry enforces degeneracy at k and 
— k + Q and allows a putative Fermi surface to be gapped 
by an infinitesimal "pairing" term with momentum Q. 
Note that a nonzero Q is encountered in the physical 
state but no gauge choice possesses one of those symme- 
tries. If the ground state has no time reversal, tt rotation, 
or inversion symmetry, a putative Fermi surface will be 
stable against any weak perturbations. In the current 
model, the tt rotation and inversion symmetries are ex- 
plicitly broken by Ja 7^ Jv or J' A 7^ J v ; the time rever- 
sal symmetry is always spontaneously broken. It follows 
that the spin Fermi surface found in our model is sta- 
ble against any weak perturbations. This is in contrast 
with the spin Fermi surface found in Refs. [TBTfTT] that 
can be gapped by some sort of weak "nesting" or "pair- 
ing" terms. Lastly, this state has recently been shown to 
exhibit bond energy correlations that have 1/|t"| 3 power 
law behaviour 37 despite having ultra-short spin corre- 
lations. 

Concluding remarks-We studied an exactly solvable 
model of spins on the kagome lattice with an odd number 
of electrons per unit cell (making it a "true" Mott insu- 
lator) that realizes a gapped chiral spin liquid with two 
gapless chiral Majorana edge modes and a gapless phase 
with a finite spin Fermi surface that is stable against 
any weak perturbations. While in terms of the spin vari- 
ables the model may seem highly anisotropic and there- 
fore somewhat artificial, Kitaev-type models have been 
derived as effective low-energy theories from more "real- 
istic" Hamiltonians with strong spin-orbit coupling [35] 
and can be readily engineered in optical lattices [331 HU] • 
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